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ABSTRACT 
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Context. Compact circumstellar envelopes (CSEs) have now been detected around several Cepheids by interferometry. These en- 
velopes are particularly interesting for two reasons: their presence could impact the Cepheid distance scale, and they could imply 
stellar mass loss. 

Aims. Our goal is to establish the spatial and spectral properties of the CSEs of t Car and RS Pup. This is done through a parametriza- 
Ph tion of the envelopes in terms of fractional flux (with respect to the star) and angular size. 

Methods. We retrieved archival Spitzer images of the two stars (A = 3.5 - 70/um), and obtained new diffraction-limited imaging with 
the VLT/VISIR camera in BURST mode (A = 8.6- 1 1.9 pm), as well as interferometric observations with VLTI/MIDI (A = 8-13 /an). 
This combination of single-telescope and interferometric techniques allows us to probe the envelopes of the two Cepheids over a broad 
range of angular scales, from arcminutes to milliarcseconds. 

Results. The circumstellar envelope of RSPup is resolved spatially at 24 and 70yimr by Spitzer, and around lOyum by MIDI and 
VISIR. The envelope of £ Car appears much more compact, and is resolved only in the VISIR and MIDI observations. The infrared 
excesses we detect around RS Pup and £ Car are both very significant, but differ considerably in spectral and spatial properties. We 
detect a warm component in the CSE of both stars at a spatial scale of a few 100 to a few 1 000 AU. In addition, RSPup presents a 
very large (several 100 000 AU) and cold (« 40 K) dusty envelope. 

Conclusions. The observed properties of the CSEs lead us to propose that the cold dust content of the large reflection nebula sur- 
rounding RS Pup has an interstellar origin, while the warm CSEs of the two stars were created by ongoing stellar mass loss. We also 
speculate that the NGC 7023 reflection nebula surrounding the Herbig Be star HD 200775 is an analogue of RS Pup at an age of 
oo : 100000 years. The presence of CSEs around the two brightest long-period Cepheids indicates that warm CSEs is probably common 

around Cepheids. However, very large dusty envelopes such as that of RS Pup are probably rare, as according to our scenario, they 
require the presence of a high dust density in the interstellar medium at the time of the formation of the Cepheid progenitor. 

Key words. Stars: circumstellar matter; Stars: variables: Cepheids; Stars: individual: £ Car, RS Pup; Techniques: high angular reso- 
lution; Techniques: interferometric; Infrared: stars 
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Ctf ■ nary of its first calibration (Leavitt & PickeringHHl, Cepheids nebula ^ but the physica i association is uncertain. The presence 

appear as increasingly complex objects. The first IRAS obser- of CSEs around Cep heids raises questions about the m ass los s 

vations in the 1980s led to the discovery of infr ared ( IR) ex- mec hanisms of these stars. As shown by Merand et al. (120071 

cess arou nd ma ny Cepheids (McAlary & Welch [1986J Deasy long period Cepheids tend to show brighter CSEs in the near-IR 

& Butler 1986|. Ba sed on these observations and IUE ultravi- K band than short per iods, indicating that the mass-loss mech- 

olet spectra, Deasy dl988J identified mass loss in a number of anism could be linked to the pu i sation of the star. Intere stingly , 

Cepheids. A very significant mass loss rate is attributed by this pu i sat i on driven mass loss models by Neilson & Lester ((2008a) 

author to RS Pup (10 6 M^yr 1 ). Bohm-Vitense & Love H1994J predict lower mass loss rates for long periods compare d to short 

obtained an even higher value of 2. 10 5 M Q .yr for /'Car. Recent per iods. Using the same models, Neilson et al. (2008) predicted 

interferometric observations revealed compact circumstellar en- a sufficiently large IR excess to impact the near and mid-IR P-L 

velopes around several nearby Cepheids (Kervella et al. 2006, relations of LMC Cepheids. 
Merand et al. 120061; Merand et al. 120071 ). Before these obser- 
vations, only one Cepheid was known with certainty to be sur- ln the present article, we report our observations of the 

two long-period Cepheids i Car (HD 84810, SAO 250683) and 

Send offprint requests to: P. Kervella RSPup (HD 68860, SAO 198944) using single-telescope and 
* Based on observations made with ESO telescopes at Paranal interferometric techniques in order to search for and character- 
Observatory, under ESO programs 073.D-0142(B), 077.D-0500(A), ize the mid-IR excess of these stars. We detail the data sets ob- 
078.D-0739(B), 078.D-0739(D) and 081.D-0165(A). tained with the Spitzer, VLT/VISIR and VLTI/MIDI instruments 
Correspondence to: Pierre.Kervella@obspm.fr in Sect. [2] Sections[3]to[5]are dedicated to the analysis of the data, 



1 . Introduction rounded by a circumstellar envelope (CSE), RS Pup. This large 

nebula (» 2') scatters the light from the central star, and hosts 
Exactly a century after the discovery of their Period-Luminosity one of the most str iking examples of light echoes known to date 
(P-L) relation (Leavitt 1908 ), and a few years from t he centen- (Kervella et al. 120081 . Another Cepheid, SU Cas, is close to a 
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while we summarize in Sect. [6] the VINCI data of I Car, previ- 
ously published by Kervella et al. (2006). We present in Sect. [7] 
the properties of the excess emission present in the spectral and 
spatial energy distribution of the two stars. The current standard 
explanation of the presence of Cepheid envelopes is that they 
were created by mass loss from the the stars themselves, through 
a mechanism probably associated with pulsation. We propose an 
alternative scenario for the formation of the CSE of RS Pup in 
Sect. [8] as well as a critical evaluation of the existing Cepheid 
mass loss models. The detailed physical modeling of the CSEs 
is beyond the scope of the present paper and is not addressed. 

In the following, we will phase the observations using ref- 
erence epochs To (expressed in Modified Julian Date) and pe- 
riods of T = 52289.916 and P = 35.55 13 days for ^Car 
(Szabados [19891 , T = 540 90.836 and P = 41.4389days for 
RS Pup (Kervella et al. [2008l >. 

2. Observations and data reduction methods 

2.1. Spitzer observations 

We retrieved the observations of ^Car and RS Pup obtained with 
the InfraRed Array Camera (IRAC; Fazio et al. I2004l > and the 
Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 
2004; Werner et al. I2004al l from the archiveQ at the Spitzer 
Science Center (SSC). The observations of £Cai with IRAC 
took place on 13 August 2006 (MJD 53 960.984, <p = 0.00), 
and with MIPS on 19 July 2006 (MJD 53 935.328, 4> = 0.28), 
while RS Pup was observed with IRAC on 29 December 2006 
(MJD 54 098.477, = 0.18), and with MIPS on 4 December 
2006 (MJD 54 073.425, = 0.58). The IRAC and MIPS wave- 
lengths complement very well our MIDI and VISIR observations 
around 10pm towards shorter and longer wavelengths. 

The IRAC images were taken in the four channels of the 
camera: 3.8, 4.5, 5.8 and 8.0 pm. We used th e MOPEX soft- 
ware suite vl8.1.5 (Makozov & Marleau 2005) to assemble mo- 
saics of the individual images, with a pixel scale of 1 .22"/pixel. 
The MIPS images were obtained in photometry mode at 24 pm 
(silicon-based array) and 70pm (germanium-based detector). 
The pipeline recipes that were applied by the archive to con- 
vert the raw data into these BCD files are described by Gordon 
et al. (120051 1 . At 70 pm, we selected the non-filtered versions of 
the BCD images. We then used MOPEX to assemble mosaics of 
the individual images, with a pixel scale of 2.45"/pixel (close to 
the original BCD images) at 24 pm and 4"/pixel at 70 pm (about 
half of the original pixel scale). 

2.2. VISIR observations 

We used the VISIR instrument (Lagage et al. 2004), installed 
at the Cassegrain focus of the Melipal telescope (UT3) of the 
VLT (Paranal, Chile). VISIR is a mid-IR imager that also pro- 
vides a slit spectrometer. As it is a ground-based instrument, its 
sensitivity is severely limited by the high thermal background 
of the atmosphere, compared for instance to Spitzer, but its re- 
solving power is ten times higher, thanks to the 8 m diameter 
of the primary mirror. VISIR is therefore very well suited for 
our programme to search for compact, spatially resolved mid-IR 
emission around I Car and RS Pup. 

However, under standard conditions at Paranal (median see- 
ing of 0.8" at 0.5 pm), the 8 m telescope is not diffraction limited 
at 10pm (seeing ss 0.4" vs. 0.3" diffraction). Instead of a pure 

1 http://ssc.spitzer.caltech.edu/ 



Airy diffraction pattern, several moving speckles and tip-tilt usu- 
ally degrade the quality of the image (see e.g. Tokovinin, Sarazin 
& Smette 2007}. To overcome this limitation, a specific mode 
of the instrument, called the BURST mode, was introduced by 
Doucet et al. (2006 2007). Its principle is to acquire very short 
exposures (At < 50 ms), in order to keep the complete integra- 
tion within a fraction of the coherence time (w 300 ms at Paranal 
in the mid-IR). The detector is therefore read very quickly, and 
the resulting images freeze the turbulence. It is subsequently 
possible to select the best images that present a single speckle 
("lucky imaging"), and are thus diffraction-limited. Our process- 
ing consists of three main steps: 1) integer-pixel recentering of 
the chop-nod background corrected frames, 2) oversampling of 
each frame by a factor of 4 using spline interpolation (resam- 
pled scale of 19 mas/pixel), 3) fine recentering with a Gaussian 
fitting procedure, 4) averaging of the best frames with respect to 
the Strehl ratio. A more detailed description of the processing is 
presented in Kervella & Domiciano de Souza (120071 1. 

The observations of t Car and RS Pup were performed dur- 
ing the nights of 23 and 24 May 2008. Table[T]lists the sequence 
of our VISIR observations of the Cepheids and their associated 
calibrator stars. They were observed immediately before and af- 
ter the Cepheids to monitor the evolution of the sky transparency 
(photometry) and thermal infrared seeing (PSF calibration). The 
observations were obtained in three filters^: PAH1, PAH2 and 
SiC covering respectively the following wavelength domains: 
8.59 + 0.42 pm, 1 1.25 ± 0.59 pm and 11.85 ± 2.34 pm. 

For an unknown reason, possibly a temporary burst of very 
bad seeing or thin cirrus clouds, observation #c suffers from a 
low photometric SNR, and was not included in the analysis. Due 
to the lower sensitivity of VISIR in the SiC filter, we could not 
recenter the individual VISIR images for observation #p. 

2.3. MIDI observations 

MIDI (Leinert et al. l2003l Ratzka et al. 120071) is the mid-infrared 
two-telescope beam combiner of the Very Large Telescope 
Interferometer (VLTI; Glindemann et al. 2004). This instrument 
is a classical Michelson interferometer working in the astronom- 
ical N band (7. 6-13. 3pm). It can combine the light from two 
VLT Unit Telescopes (8.2 m, UTs) or two Auxiliary Telescopes 
(1.8 m, ATs). For the reported observations of t Car and RS Pup, 
we used a prism with a spectral resolution of R = A/AA ^ 30 
to obtain spectrally dispersed fringes. During the observations, 
the secondary mirrors of the two Unit Telescopes were chopping 
with a frequency of 2 Hz to properly sample the sky background. 
MIDI has two photometric calibration modes: SCI_PH0T, in 
which the photometry is measured simultaneously with the inter- 
ferences fringes, and HIGH_SENS, in which the flux is measured 
separately before the interferometric observations. The visibil- 
ity is derived from the ratio of correlated flux (modulation in 
the interferometric fringes) with the total flux (given by the pho- 
tometry). Due to the relatively low thermal-IR brightness of the 
two Cepheids, the observations presented here were obtained in 
HIGH_SENS mode, using the low dispersion setup (spectral reso- 
lution of A/AA « 30). 

For the raw data processing, we used two different soft- 
ware packages: MIA developed at the Max-Planck-Institut fiir 
Astronomie and EWS developed at the Leiden Observatory 
(MIA+EW5fl version 1.5.2) in order to extract first the instrumen- 
tal squared coherence factors and then the calibrated visibilities 

2 http://www.eso.org/sci/facilities/paranal/instruments/visir/inst/ 

3 http://www.strw.leidenuniv.nl/~nevec/MIDI/index.html 
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Table 1. Log of the VISIR observations of ^Car and RSPup. 
MJD* is the modified Julian date of the start the exposures on 
the target minus 54000, and <p the phase of the Cepheid. The 
Detector Integration Time (DIT) is given in milliseconds for one 
BURST image. 9 is the seeing in the visible {A = 0.5 jjm) mea- 
sured by the observatory DIMM sensor, in arcseconds, and AM 
is the airmass of the observation. The duration of each observa- 
tion was 360 s on target. 



# 


MJD* 




Star 


Filter 


DIT 


e 


AM 


a 


609.953 




HD76110 


PAH1 


16 


1.5 


1.09 


b 


609.961 




HD 76110 


PAH2 


8 


1.4 


1.11 


c !:! 


609.980 


0.53 


RS Pup 


PAH1 


16 


1.6 


1.26 


d 


609.987 


0.53 


RS Pup 


PAH2 


8 


1.3 


1.31 


e 


610.002 


0.26 


^Car 


PAH1 


16 


1.5 


1.34 


f 


610.009 


0.26 


^ Cai- 


PAH2 


8 


1.3 


1.37 




610.035 




ro 89682 


PAH1 


16 


1.6 


1.26 


h 


610.042 




HD 89682 


PAH2 


8 


1.8 


1.29 


i 


610.081 




HD 98118 


PAH1 


16 


1.7 


1.31 


j 


610.089 




HD 98118 


PAH2 


8 


1.6 


1.36 


k 


610.960 




HD 76110 


PAH1 


16 


0.7 


1.11 


1 


610.967 




HD76110 


SIC 


20 


0.9 


1.13 


m 


610.979 


0.29 


^Car 


PAH1 


16 


1.0 


1.31 


n 


610.986 


0.29 


^Car 


SIC 


20 


0.8 


1.32 





610.998 


0.55 


RS Pup 


PAH1 


16 


0.8 


1.38 


P 


611.005 


0.55 


RS Pup 


SIC 


20 


0.9 


1.44 


q 


611.016 




HD 89682 


PAH1 


16 


0.6 


1.23 


r 


611.024 




HD 89682 


SIC 


20 


0.9 


1.24 


s 


611.073 




HD 124294 


PAH1 


16 


0.8 


1.06 


t 


611.081 




HD 124294 


SIC 


20 


0.8 


1.05 



* Observations #c and #p suffer from low photometric SNR. 



Table 2. Log of the MIDI observations of £Car and RSPup. 
<p is the phase of the observation for the Cepheids, relative to 
the maximum light in V (see Sect. Q] for P and Tq values), B 
is the projected baseline length, PA is the position angle of the 
projected baseline (counted East of North = 0°), and AM is the 
airmass. 



# 


MJD 


<P 


Target 


B(m) 


PA(°) 


AM 


A* 


53104.060 


0.90 


e Car 


40.82 


48.7 


1.27 


B 


53104.076 




HD 67582 


38.86 


58.9 


1.20 


( 1 




n on 


c L,ar 


jy.Lo 


JO. 4 


1 


u 


J j 1U4. 1U4 




1 1U / JoZ 


JJ.OJ 


OD.J 


1 1 1 


E 


53 842.082 


0.66 


e Car 


117.52 


74.4 


1.30 


F 


53 842.099 




HD 67582 


122.80 


70.1 


1.21 


G" 


54 045.290 


0.90 


RS Pup 


121.25 


31.1 


1.29 


H 


54 046.275 


0.92 


RS Pup 


119.80 


27.5 


1.36 


I 


54 046.312 




HD 70555 


122.87 


37.2 


1.20 


j,« 


54 046.343 


0.41 


£Cm 


130.10 


26.5 


1.52 


K 


54 046.372 




HD94510 


130.11 


20.7 


1.57 


L 


54 047.267 


0.95 


RS Pup 


119.10 


25.7 


1.40 


M 


54 047.286 




HD 70555 


119.62 


30.3 


1.31 


N 


54 048.244 




HD 70555 


114.88 


17.3 


1.63 


O 


54 048.270 


0.97 


RS Pup 


119.79 


27.5 


1.36 


P 


54 048.347 


0.46 


^Car 


130.03 


29.3 


1.48 


Q 


54 048.366 




HD94510 


130.10 


20.4 


1.57 



* ( Car observations obtained in 2004 (see Sect. |2.3| l. 
** Observation #G of RS Pup was not calibrated. 
*** Observation #J of C Car has a low interferometric flux. 



V(A) (Chesneau 2007). We found a good agreement between the 
results of the MIA and EWS packages within the error bars. In 
the following we will refer to the results obtained with the EWS 
package. To remove the instrumental and atmospheric signa- 
tures, three calibrators were observed together with the Cepheids 
: HD 70555 (G LD = 2.54 + 0.04 mas) for RSPup, HD 67582 
(6 LD = 2.39 + 0.06 mas) and HD 94510 (0 LD = 2.23 + 0.0 1 mas) 
for ^Car. They were chosen in the Cohen et al. (1999) cata- 
log, except HD 94510 which was selected from the CalVirQ 
database. The calibrators are nearly unresolved by MIDI on the 
selected baselines (V > 95%), and the systematic uncertainty 
associated with their a priori angular diameter error bars is neg- 
ligible compared to the typical precision of the MIDI visibilities. 

The log of the MIDI observations of £ Car and RS Pup is pre- 
sented in Table [2] The new observations took place on 17 April 
2006, and 7-9 November 2006, on the UT1-UT4 baseline. To 
complete this new data set, we retrieved from the ESO archive 
the observations of I Car obtained on the night of 8-9 April 2004 
on the UT2-UT3 baseline, and we reprocessed them using the 
EWS package. These observations were originally presented by 
Kervella et al. (2006). but they were processed at the time us- 
ing the MIDI Data Reduction Software developed by the Paris 
Observatory. The present reprocessing is intended to ensure the 
consistency of the visibilities with the more recent observations. 
The 2004 observations were obtained with natural seeing, while 
the 2006 observations employed the MACAO adaptive optics 
system (Arsenault et al. 2003), that provides an improved sta- 
bility of the photometry in MIDI. 



4 http://www.eso.org/observing/etc/ 



The mean photospheric angular diameters of /'Car and 
RSPup are respectively #ld = 2.992 + 0.012 mas (measured 
by Kervella et. al. 2004ajl, and 6>ld ~ 1.015 mas (predicted by 
Moskalik & Gorynya 120051 |2006b . Considering the projected 
baseline lengths and the wavelength of the MIDI observations, 
the photosphere of /'Car is only slightly resolved at 8 /jm with 
V ~ 93% and RSPup is mostly unresolved, with V ~ 99%. 
A measured visibility significantly lower than these figures will 
therefore be characteristic of the presence of additional emission 
close to the star. The analysis of the MIDI visibilities is presented 
in Sect. 15 ,4| and |5 .5 [ respectively for fCai and RS Pup. The pres- 
ence of a strong ozone atmospheric absorption band over the 
range A = 9.3 - 10.0;um (Lord 1992) can make the visibilities 
unreliable in this domain, as the photometric calibration is made 
more difficult by the lower flux level. 

2.4. Fourier analysis technique 

In Sections [3] |4]and|5] we present the analysis of the spatially re- 
solved emission around I Car and RS Pup as observed by Spitzer, 
VISIR and MIDI. The spatial scale sampled by these three in- 
struments covers a broad range from arcminutes with Spitzer 
to arcseconds with VISIR and milliarcseconds with MIDI. This 
coverage allows us to probe the different scales of the infrared 
emission around the two Cepheids, both angularly and spec- 
trally. 

In the present work, we search for spatially extended emis- 
sion in the Spitzer and VISIR images using a Fourier based tech- 
nique, similar in its principle to the calibration technique used 
for interferometry. For this purpose, we divide the Fourier trans- 
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form modulus of the image of the Cepheid (/cepheid) by that of 
the calibrator image (/calibrator): 



¥(v*,v y ) 



FFT[/ C epheid(^,3 ; )] 



FFT [/, 



Calibrator 



ix,y)] 



(1) 



where (x,y) and (v x ,v y ) are conjuguate Fourier variables (sky 
coordinates and angular spatial frequencies for instance). This 
Fourier technique has three major advantages over direct image 
analysis techniques (e.g. PSF subtraction). Firstly, it is similar 
in its principle to a deconvolution of the image, but the absence 
of an inverse Fourier transform step after the division in Eq. Q] 
avoids the appearance of uncontrolled phase noise. Secondly, 
the background noise is rejected more efficiently, as it is well 
separated from the image itself by its spatial frequency spec- 
trum. It is therefore more sensitive than direct image analysis 
techniques. Thirdly, interferometric observations provide mea- 
surements of the Fourier transform of the intensity distribution 
of the observed object. The Fourier analysis applied to images 
allows us to analyze them within exactly the same framework as 
interferometric observations. The derived CSE properties (an- 
gular size, flux contribution) are therefore strictly comparable, 
which is particularly interesting for the present multi-instrument 
study. 

We compute the ring median of ^ (i.e. median for a given 
spatial frequency radius v over all azimuth directions in the 
Fourier space) to increase the SNR. This operation provides 
us with a quantity *P(v) that is mathematically equivalent to 
the squared visibility of interferometric fringes. In the present 
Fourier-based analysis, we will always consider ring median 
properties in order to obtain the global properties of the CSEs. 
The zero frequency, corresponding to the average ratio of the 
Cepheid and calibrator images, is normalized to unity in all 
cases. Our estimation of the error bars on *P takes into account 
both calibration and statistical uncertainties. The calibration un- 
certainties are derived from the relative dispersion of the PSF 
calibrator's Fourier modulus over the night. The statistical uncer- 
tainty is taken as the rms dispersion of the calibrated ^ function 
over all azimuth directions, for each spatial frequency. Both error 
sources are added quadratically to get the total uncertainty. Our 
choice to consider azimuth-averaged values implies that devi- 
ations from central symmetry of the CSEs will be undetectable. 
However, the error bars associated with the derived ¥ values 
include the dispersion that may arise from asymmetric flux con- 
tributors. 

From the Van Cittert-Zernike theorem, the *P(v) functions 
are the Fourier transforms of the spatial flux distribution of the 
CSEs. It is therefore possible to retrieve the CSE intensity distri- 
bution from the shape of these curves, by adjusting the Fourier 
transform of an approximate model of the intensity distribution. 
We thus define the simple model of a point-like star surrounded 
by a Gaussian shaped CSE with a FWHM of p,\ and a relative 
flux contribution of 



/cse(^) 
/*(<*) ' 



(2) 



where /cse is the photometric flux from the Gaussian envelope 
and f* the flux from the unresolved component. This model is 
suitable for the Spitzer, VISIR and MIDI UT2-UT3 observa- 
tions, for which the FWHM of the stellar photosphere is sig- 
nificantly smaller than the diffraction limit A/D of the instru- 
ment (with D the aperture). For the longer baseline UT1-UT4 
observations with MIDI, the photospheric extension must be 



90% - 

&■ 

CO 

(3 80% - 




• Star+0% CSE 
a Star+10% CSE 1" 
Star+30% CSE 1 " 

Model 0% 

Model 10% 

— Model 30% 
- — Unresolved 



0,0 0,5 1,0 

Spatial frequency v (1/arcsec' 

Fig. 1. *P function in the PAH1 filter resulting from the Fourier 
analysis of HD 89682 considered as the scientific target, with 
HD98118 and HD 124294 as PSF calibrators. To validate our 
Fourier analysis method, we artificially introduced two Gaussian 
CSEs with fluxes of 10 and 30% of the stellar flux and a FWHM 
of 1", as discussed in Sect. 12.41 

taken into account (see Sect. |5]l. This type of model was already 
used by Kervella et al. (2006) and Kervella & Domiciano de 
Souza (2006). The l'(v) function of this model is: 



1 + a A 



1 + at exp 



(npvf 
' 4 In 2 



(3) 



2.5. Fourier technique validation 

As a validation test for our Fourier analysis technique, we ap- 
plied it to the VISIR observations of HD 89682 obtained in the 
PAH1 filter, using it as the scientific target (observations #g and 
#q in TableQ). We used HD 981 18 and HD 124294 as PSF cal- 
ibrators (observations #i and #s). As shown in Fig. Q] (upper 
curve), the resulting average function of HD 89682 is sta- 
tistically identical to 100%, as expected for an unresolved tar- 
get. The derived flux contribution for the extended emission is 
a = 0.4 + 0.5%, while the value of p is not constrained. 

To check the sensitivity of the Fourier technique to the pres- 
ence of a resolved Gaussian CSE component, we created ar- 
tificial images using a point-like central source surrounded by 
an extended Gaussian emission with a fixed FWHM of 1" and 
two values of its flux contribution a relative to the central star: 
10% and 30%. These model images were then convolved with 
the #g and #i images of HD 89682 separately, to obtain real- 
istic images (separately for each observing night) of the model 
as observed with VISIR. We then processed these synthetic im- 
ages separately for each night using the Fourier technique. We 
finally took the average of the functions of the two nights, that 
are shown in Fig. Q] (lower two curves). By fitting the Gaussian 
model presented in Eq. [3] we finally derive the following a and 
p parameters for the star+10% CSE model: 



p(Model 10%) = 0.99 + 0.12 arcsec, 
a(Modell0%) = 11.1 ± 1.0%, 
and for the star+30% CSE model: 
p(Model 30%) = 1.00 ± 0.04 arcsec, 



(4) 
(5) 



(6) 
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Table 3. Summary of the results of the analysis of the Spitzer 
images. 



Technique 


IRAC 


MIPS 24 pm 


MIPS 10 pm 


£Car 


Image 


saturated 


unresolved 


unresolved 


Fourier 


saturated 


unresolved 


unresolved 


RSPup 


Image 


unresolved 


p * 26" 
a « 100% 


p * 118" 
(star invisible) 


Fourier 


unresolved 


p = 60 ± 3" 

a- = 98.6± 1.3% 


not applicable 
(star invisible) 



a(Model30%) = 31.2 + 1.2' 



(7) 



The derived relative fluxes and FWHMs are in agreement within 
lcr with the input model values, thus validating our Fourier- 
based approach. In the following Section, we detail the appli- 
cation of this analysis approach to the Spitzer, VISIR and MIDI 
instruments, by order of increasing resolution. 




3. Spitzer data analysis 

In this Section, we present several analysis methods applied to 
the Spitzer images of t Car and RS Pup. Table [3] summarizes the 
results of these different techniques. 




Fig. 2. Spitzer images of the reference PRF and £Car at 24 and 
70 Ami. 



3.1. { Car images 

I Car was observed with both the IRAC and MIPS instruments 
of Spitzer, in all the available spectral channels of these instru- 
ments: 3.8 to 8.0yum for IRAC, and 24 to 160 pm for MIPS. 
The 160/^m observation is not useable, and the IRAC images 
of /'Car are all saturated close to the star position, and therefore 
are not suitable for small-scale CSE searches. On a larger spa- 
tial scale, the IRAC images do not present an obvious extended 
emission that would indicate the presence of a very extended 
nebula (> 30" from the star). 

In order to check for the presence of spatially extended emis- 
sion from the immediate surroundings of (Cai in the MIPS im- 
ages, we used as references the experimental point response 
functions (PRFs) provided in the MOPEX calibration pack- 
age that are elaborated from calibration observations of the 
Trapezium in Orion. We also checked synthetic PRFs from the 
sTinyTim software packag^fl but they do not provide as clean 
subtraction residuals as the MOPEX experimental PRFs. 

The MIPS images of the PRF reference and I Car at 24 and 
70/iin are shown in Fig. [2] The two images on the left show the 
oversampled PRF of Spitzer (0.6125"/pix) and the 24 pm image 
of ^Car (2.45"/pix). The two images on the right show the same 
data for 70 pm (pixel scales of 1 and 4"/pix). The square-root 
color scale at 24 pm has been chosen to enhance the faint exten- 
sions of the PRF, hence it cuts the inner core of the star images. 
The 70 pm color scale is linear, with contours corresponding to 
3 to 30% of the peak level (steps of 3%). 

3.2. RS Pup images 

RS Pup was observed both with IRAC's sub-array and MIPS (in 
all channels except at 160 pm). We searched for extended emis- 
sion in the IRAC images using the Fourier transform technique 

5 http://ssc.spitzer.caltech.edu/archanaly/contributed/stinytim/ 



presented in Sect. 12.41 We used for this purpose a PSF calibrator, 
BD+68 1022, that was observed on 20 June 2008 specifically 
for the calibration of the instrumental PRF. RS Pup appears un- 
resolved in the four channels of IRAC. From the resolution of the 
IRAC observations, we can set an upper limit on the extension 
of the CSE to 5" in the IRAC bands. 

The left part of Figure [3] shows the circumstellar nebula of 
RS Pup in the V band, from observations with the NTT/EMMI 
instrument by Kervella et al. (2008). In this wavelength range, 
we observe the scattered light from the central star. The middle 
and right images are the Spitzer/MlPS images at 24 and 70 pm. 
They show the thermal emission from the circumstellar dust (the 
three images have the same spatial scale and orientation). The 
nebula appears elongated along a north-east axis both in scat- 
tered light and IR emission. Interestingly, RS Pup itself is offset 
from the brightest part of the nebula, particularly at 70 pm. It 
also remains invisible at this long wavelength as the image is 
fully dominated by the thermal emission from the envelope. 

3.3. Radial CSE intensity distribution 

The average surface brightness at 24 and 70 pm of RS Pup's CSE 
as a function of the distance from the central star is shown in 
Fig. |4] The maximum surface brightness of the CSE of RS Pup 
reaches 0.1 mJy.arcsec" 2 at 24 pm (outside the star image) and 
more than 4 mJy.arcsec -2 at 70 pm. While the CSE is inhomo- 
geneous locally, its overall surface brightness follows closely an 
exponential law as a function of the angular separation from the 
star. An ad hoc fit from 20 to 90 " from the star at 24 pm gives: 

/24„ m (RSPupCSE) = 201.6 exp(-^) (8) 

with f24/ lm in /uJy/arcsec 2 and 6 in arcsec. This corresponds to a 
full width at half maximum of 26 " at 24 pm (but the distribution 
is not strictly Gaussian). As the wings of the PSF of RS Pup 
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Fig. 3. Images of the nebula surrounding RS Pup in the V band (left, pulsation phase (p = 0.877), at 24 //m (middle) and at 70 fim 
(right), to the same spatial scale. The images are not corrected for the background emission. The position of the star is marked with 
the * symbol. The contours are 576 to 705 /dy.arcsec -2 (step 5.9) in the 24 /mi image, and 627 to 5012/dy.arcsec~ 2 (step 627) in 
the 70/mi image. The linear scale is given for a distance of 1 992 pc (Kervella et al. 120081 1. 



10000,0 



1000,0 




0,1 J , , , , , , , , , , , , 

10 20 30 40 50 60 70 80 90 100 110 120 130 
Angular radius from RS Pup (arcsec) 

Fig. 4. Surface brightness of the CSE of RS Pup at 70 /mi (top, 
black curve) and 24 fmx (bottom, grey curve) as a function of the 
angular radius to the star. The result of an exponential fit to the 
24 /mi curve between 20 and 90 arcsec is shown as a dotted line. 
The dot-dashed line is a Gaussian fit to the 70/im profile (see 
text for details). 



3.4. Fourier analysis 

An alternative method to analyse the IRAC and MIPS images is 
through a Fourier analysis, using the model defined in Sect. I2.4I 
For the calibration of the *P function of the RS Pup images ob- 
tained with IRAC, we used the archive images of BD+68 1022 
that was observed on 20 June 2008 specifically for the cali- 
bration of the instrumental PRF. The resulting functions for 
the four IRAC channel images show that RS Pup appears unre- 
solved. We cannot apply the Fourier analysis technique to the 
t Car images as they are saturated close to the star. 

We analyzed the MIPS images using the theoretical PRF ref- 
erence presented in Fig. [2] The result of the T* function cali- 
bration is presented in Fig. [5] ^Car appears unresolved in the 
MIPS images at 24 and 70 //m, with a *P function close to 100% 
at all spatial frequencies. As expected from the visual aspect of 
RS Pup in the 24 /mi MIPS band (Fig.|5), the *F function shows 
a steep decline around a spatial frequency of v = 0.015 arcsec -1 . 
This frequency corresponds to the typical w 1' angular radius of 
the nebula, as shown in Fig. [3] In this case, the adjustment of our 
Gaussian CSE model using a classical^- 2 minimization gives: 

P24^m(RS Pup) = 60 ± 3 arcsec, (10) 
ff 24pm (RS Pup) = 98.6+ 1.3%. (11) 



are included in this fit, they tend to create a more peaked light 
distribution than the nebular emission alone. 

The distribution of flux in the nebula at 70 /mi does not show 
this effect, as the stellar PSF is negligible compared to the enve- 
lope emission. It closely follows a Gaussian law of the form: 



/tgWRS Pup CSE) = 1081 exp 



-I—)' 
V49.6/ 



(9) 



with fjopm in /dy/arcsec 2 and 9 in arcsec. This corresponds to a 
full width at half maximum of 1 18 ". 



We here retrieve a larger FWHM than from the direct exponen- 
tial model-fitting to the image shown in Fig. [4] This is due to the 
fact that the Fourier technique cleanly separates the unresolved 
stellar emission, including the wings of the stellar PSF, and the 
resolved nebular emission. 

The analysis of the 70 /mi MIPS image of the nebula is better 
done directly on the images using a simple Gaussian fit (as done 
in Sect. I3.3K because the star is invisible and the nebula is well- 
resolved in the image. 

3.5. Aperture photometry 

We explored two routes to obtain the photometry of i Car and 
RS Pup from Spitzer images: using the APEX single-frame tool 
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Fig. 5. Calibrated *F function of the MIPS images of £ Car and 
RS Pup. 



included in MOPEX to obtain PRF-fitting fluxes, and from clas- 
sical aperture photometry using self-developed routines. The re- 
sults of the two approaches are in excellent agreement for £ Car, 
but it appeared that due to the presence of the irregularly shaped 
CSE around RS Pup in the MIPS images, the PRF-fitting was not 
robust on this star. For this reason, we adopted the aperture pho- 
tometry technique for both stars at MIPS wavelengths. For IRAC 
images, we adopted a slightly different measurement technique 
for /'Car and RS Pup due to the saturation of the core of £ Car's 
images (discussed below). The PRF-fitting method is in any case 
not recommended by the Spitzer Science Center for IRAC im- 
ages. The systematic uncertainty in the conversion of the IRAC 
instrumental fluxes to absolutely calibrated units is estimated by 
Reach et al. (2005) to be less than 3% (value retained here). 
The same values for MIPS fluxes are found by Engelbracht et 
al. (I2007I) to be 4% at 24 pm, and by Gordon et al. d2007l > to be 
5% at 70/im. As we observe bright stars, we neglect the statis- 
tical uncertainties on our photometry and adopt these values as 
the total uncertainties of our measurements. The measured flux 
densities on the two stars are summarized in Table [5] while the 
details of the photometry are given in the following two para- 
graphs. 

3.5.1. IRAC 

Due to the heavy saturation of the core of the IRAC images of 
£ Car, we chose to compare the fluxes of the theoretical extended 
PRF referenc^ and /'Car over an unsaturated ring between an- 
gular distances of 10 and 20 pixels (1.22 "/pixel) from the star. 
This angular separation range allows us to stay in the linear part 
of the detector's sensitivity curve: we have a sufficient SNR for 
the measurement, but avoid a bias from the saturated core. For 
an example of core -saturated images obtained with IRAC, the 
reader is referred to Marengo et al. (2006). We corrected the 10- 
20 pixel flux of £ Car for the sky background using an annulus of 
20-30pixels, in the four IRAC channels. We measured the rele- 
vant aperture corrections y on the extended PRF image for this 
specific combination, according to the IRAC Data Handbook 



6 http://ssc.spitzer.caltech.edu/irac/psf.html 



version 3.C0 We then obtained the total flux of the star from 
a multiplication by y. 

For RS Pup, the photometric analysis is easier than at longer 
wavelengths as the star appears unresolved (Sect. 13. 2t . and we 
obtained classical aperture photometry. Due to the small field of 
view, we used a small aperture of 5 pixels, and a background 
annulus between 5 and 10 pixels. The aperture corrections listed 
by the SSC for these apertures are 1.061, 1.064, 1.067 and 1.089 
respectively for the 3.8, 4.5, 5.8 and 8.0/mi channels. 

3.5.2. MIPS 

The aperture photometry of fCar is simplified by the fact that it 
appears as a bright (but unsaturated) isolated point source within 
a relatively large field void of other bright sources. Taking ad- 
vantage of this configuration, we obtained the aperture photom- 
etry of ^Car at 24 pm in a classical way, using a 35" aperture 
radius, and a sky background annulus of 40-50". This configu- 
ration corresponds to an aperture correction of 1 .082 according 
to the SSC0. For the 70 pm image, we used a 30" aperture, with 
a sky annulus of 40-60", for an aperture correction of 1.295. 

The CSE of RS Pup has a significant extended contribution 
at 24pm and it fully dominates the 10pm Spitzer image. The 
goal of the present work being to establish the total infrared ex- 
cess of the star, the circumstellar and photospheric emissions are 
considered together in the photometry. For this reason, we used 
a large photometric aperture of 3.0' in radius at 24 pm to include 
the full CSE emission. At 70 pm, a difficulty is that the nebula 
appears larger than the field of view of the MIPS mosaic along 
its narrower direction, as shown in Fig. [3] As a consequence, 
there is a slight truncation of the CSE faint extensions. For this 
reason we chose to integrate the average Gaussian profile given 
in Eq. [9] from zero to infinite radius to derive the total nebular 
flux density, and obtain /7o^m(RS Pup CSE) = 16.9 Jy. 

We can estimate the compact emission at 24 pm by extrapo- 
lating the flux over a small aperture that includes essentially the 
star and little of the extended nebula. Using an aperture of 13" 
and an aperture correction of 1.167, we obtain a flux of 0.39 Jy 
for the core of RS Pup, or about half of the total flux of the star 
and its CSE. An analysis of the radial flux distribution in the 
nebula is presented in Sect. 13.41 



4. VISIR data analysis 

In this Section, we present several analysis methods applied to 
the VISIR observations of £ Car and RS Pup. Table|4]summarizes 
the corresponding results. 

4. 1 . Cepheid and calibrator images 

The VISIR images of £Car in the PAH1, PAH2 and SiC filters 
are shown in Figs. [6] and [7] together with their median radial pro- 
files. These 1(6) profiles were obtained by computing the median 
of the pixels located within an annulus with internal radius 6 and 
external radius 0+1 pixel (resampled scale of 19 mas/pixel). 

The images of RS Pup in the PAH1 and PAH2 bands are pre- 
sented in Fig. |8] The observation #n (Table [TJ of RS Pup in the 
SiC filter did not give a satisfactory image due to the limited 
sensitivity of the instrument in this band. 



7 http://ssc.spitzer.caltech.edu/irac/dh/ 

8 http://ssc.spitzer.caltech.edu/mips/apercorr/ 
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Fig. 6. VISIR #i image of {Car in the PAH1 filter (left) and median radial profiles (normalized to unity at maximum) for the two 
observations #i and #m (see Tabled, with the associated PSF references. 



Table 4. Summary of the results of the analysis of the VISIR 
images. 



Technique 


PAH1 


PAH2 


SiC 


{Car 


Image 


resolved 


resolved 


unresolved 


PSF sub. 


no residual 


no residual 


unresolved 


Fourier 


p = 0.75+0.15" 
a = 4.0 ± 1.0% 


p = 1.17 ±0.18" 
a = 10.6 ± 1.2% 


unresolved 


RSPup 


Image 


resolved 


resolved 


too faint 


PSF sub. 


faint residual 


residual x 10% 


too faint 


Fourier 


p = 0.8 + o j arcsec 
a = 2.3 ± 1.2% 


p= 1.04 ±0.21" 
a = 23.5 ±2.0% 


too faint 



4.2. PSF subtracted images 

The RS Pup image in the PAH2 band shows some asymmetry. 
An additional emission appears in the south-western quadrant 
relative to the star. In order to remove the contribution from the 
star itself, we subtracted the average image of the PSF calibra- 
tors, scaled using a classical x 2 minimization of the residuals. 
The resulting image is shown in Fig. [9] (right). This emission 
is also apparent in the original image (Fig. |S), specifically as 
an asymmetry of the first Airy ring that appears brighter in its 
south-western section. The integrated flux density of this addi- 
tional component is 2.8 10 -15 W/m 2 /yum, or « 10% of the total 
emission of RS Pup in this band. Its position is 0.4"from the 
star, along an azimuth of 220° and its FWHM angular extension 
is approximately 0.5". 

Although it appears slightly more extended than a point 
source (for which the FWHM would be 0.3"), it is not excluded 
that this additional component is stellar in nature. Its position 
on the first Airy ring of RS Pup could have distorted its shape. 
However, it is more likely a particularly dense clump of dust 
in the close-in CSE of RS Pup. Interestingly, its position places 
it almost aligned with the major axis of the nebula at 10 pm as 
shown in Fig. [3] 

The PSF-subtracted PAH1 image (Fig.|£l left) also shows a 
small extended emission in the same quadrant as in PAH2, but its 
detection is more uncertain. A compact flux contributor is also 
present to the East of RS Pup, but its detection is very uncertain 
due the high PSF residual noise level so close to the star. 



PAH2- 11.25 Mm 




Fig. 7. VISIR #f and #n images of {Car in the PAH2 and SiC 
filters (normalized to unity at maximum) and the corresponding 
median radial profiles (with the PSF calibrators). 



We applied the same processing to the images of i Car, but 
they do not reveal any significant subtraction residual. 



4.3. Fourier analysis 

We search here for extended emission in the VISIR images us- 
ing the same Fourier transform technique as for the Spitzer data. 
A spatially extended emission appears as a deficit in the *F func- 
tion, starting at the spatial frequency corresponding to the angu- 
lar scale of the emission. As for Spitzer, we consider here ring 
median values of *F(v) over all azimuth directions. It is impor- 
tant to note that the derived a contributions correspond to the 
extended emission contained within the full field of view of the 
VISIR images (3 .2" x 3 .2") and extrapolated outside this field as 
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Fig. 10. Fit of a Gaussian CSE + unresolved point source model to the *P functions of ^Car and RSPup in the VISIR PAH1 
(A = 8.59 /mi, left) and PAH2 (A =11.25 yum, right) bands. 




Fig. 8. VISIR image #o and #d of RS Pup in the PAH1 and PAH2 
filters (left) and the corresponding median radial profiles com- 
pared to the PSF calibrators. 

a Gaussian. In the PAH1 band (Fig. [10] left), t Car and RS Pup 
present a similar level of resolution, with: 



P$.59iim({ Car) = 0.75 + 0.15 arcsec, 
a8.59/im(£Car) = 4.0 + 1.0 %, 
P8.59 A ,m(RS Pup) = 0.8!°;* arcsec, 
ar 8 .59pm(RSPup) =2.3 ±1.2%. 



(12) 
(13) 
(14) 
(15) 



In the PAH2 band (Fig. [TO] right), the *P functions also show 
that the two stars present a spatially resolved flux contributor. 
The best-fit parameters for I Car are: 



RS Pup - PSF Subtracted - 8.59 nm 

N| 

1 arcsec E | 


RS Pup - PSF Subtracted - 11 .25 

N| 

1 arcsec E | 







Fig. 9. PSF subtracted VISIR images of RSPup in the PAH1 
(left) and PAH2 (right) bands. The position of the star is marked 
with 



arii.25pn(^Car) = 10.6 ±1.2%, 



(17) 



The *P function of RS Pup in the PAH2 filter shows a "bump" 
for spatial frequencies around 0.5". This behavior could be due 
to the compact emission source present in the field of view of 
VISIR (Sect. l4~2l . The Gaussian model fitting gives the follow- 
ing properties for RS Pup: 



Pn.25 /J m(RS Pup) = 1.04 + 0.21 arcsec, 
aii.2 5/I m(RS Pup) = 23.5 ± 2.0 %. 



(18) 
(19) 



Pn.25//m(^Car) = 1.17 + 0.18 arcsec, 



(16) 



The SiC image of £ Car does not show a resolved compo- 
nent. The broad filter bandpass (40%) degrades the contrast of 
the Airy diffraction pattern, as different wavelengths correspond- 
ing to different spatial resolutions are superimposed. As the de- 
tected CSE in PAH1 and PAH2 is compact, this likely prevents 
its detection. 



4.4. Aperture photometry 

We derive here the photometric flux density of t Car and RS Pup 
in the VISIR filters through classical aperture photometry. In or- 
der to absolutely calibrate the flux density, we used as references 
the spectrophotometric templates from Cohen et al. ( 1999). The 
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irradiance of the calibrators was read from the templates tak- 
ing into account the transmission curve of each filteo The in- 
tegrated fluxes are affected by differential atmospheric absorp- 
tion due to the different airmass of each observation (listed in 
Table. [TJ. To take this effect into account, we used the empirical 
formula by Schiitz & Sterzik d20051 l. that gives the multiplica- 
tive correction C(A, AM) to apply to the raw flux to calibrate the 
atmospheric absorption: 



C(A, AM) = 1 + 



0.104 

0.220 0i-8.6jum) 



(AM-1). (20) 



As we used the same apertures for the photometry of the 
Cepheids and the flux calibrators (1.3" in angular radius), no 
aperture correction is required. 

The measured flux densities on t Car and RS Pup in the 
VISIR filters are summarized in Table [5] The angular radius of 
the effective aperture over which the VISIR photometry was ob- 
tained is 1.31" in all filters. The associated uncertainties include 
the statistical dispersion of the aperture photometry, the disper- 
sion of the calibrator flux densities over the night, and the abso- 
lute calibration uncertainty (^ 1.8% for our calibrators) as listed 
by Cohen et al. (fl999l >. 

5. MIDI data analysis 

5. 1 . Differential visibilities 

The absolute calibration of the MIDI visibilities has a typical 
accuracy of 10-15% in the HIGH.SENS mode. The SCI_PH0T 
mode allows for a better accuracy, but is only available for very 
bright objects. As our two Cepheids are relatively faint stars 
for MIDI standards, the HIGH_SENS mode was the only avail- 
able choice. Its large calibration uncertainty would completely 
mask the small contributions of the CSE in the absolute visibil- 
ity spectrum. However, the shape of the spectrum is known with 
a much better accuracy, and the relative values of the visibil- 
ity between different wavelengths are accurate to a few percent 
(Chesneau 2007). In the following, we therefore adopt a differ- 
ential approach for the interpretation of the MIDI visibilities. We 
define the differential visibility as: 



AV(A,B) = 



V(A, B) 
V(8 pm, B) 



- 1 



(21) 



with B the projected baseline of the measurement, and A the 
wavelength. By definition, AV(8 pm, B) = 0. In practice, we es- 
timated the value of V(8 pm, B) from a linear fit to the whole 
visibility spectrum. To estimate the uncertainties, we considered 
the standard deviation of the residuals of this fit to AV, uniform 
with wavelength. This differential visibility approach gives us a 
very much improved sensitivity to color-dependent variations of 
the CSE properties. However, it comes with a price, as the sig- 
nature of any achromatic component is obliterated (Sect. 15.31 1. 
In the following, we compare the resolved flux component mea- 
sured with VISIR in the PAH1 and PAH2 bands to the MIDI dif- 
ferential visibilities, taking into account the expected differential 
visibility contribution from the stellar photosphere, as discussed 
in the next paragraph. 

5.2. Visibility model 

The MIDI differential visibility measurements AV{A) character- 
ize physically the ratio of the spatially correlated flux / CO n (^, B) 



(measured from the contrast of the interference fringes) to the 
total flux f iot (,A), at the spatial frequency A/B sampled by the 
interferometer, relative to this ratio at 8 pm. The total flux con- 
sidered here is the flux enclosed within the aperture stop of the 
instrument. For MIDI , this aperture has an angular radius of 
approximately 0.4arcsec (see also Sect. I5.61 >. Within this aper- 
ture, the flux from the resolved component of the CSE around 
the Cepheids contributes to increase / tot up to the edge of the 
aperture stop, even beyond the core of the point spread func- 
tion (PSF) of the telescopes. This is different for instance from 
fibered interferometric beam combiners (e.g. VINCI; Kervella et 
al. 1200 4b), whose field of view closely matches the PSF. 

The goal of this analysis is to search for a compact compo- 
nent in the CSE, in addition to the extended component resolved 
with VISIR. Our approach is the following: we present here- 
after a model of the known differential visibility contributors, 
and compare its prediction with the observed values. In this pro- 
cess, we will initially suppose that the CSE flux density at 8 pm 
is negligible, i.e. a(8 pm) = 0. This hypothesis is reasonable in 
view of the small photometric infrared excess measured at this 
wavelength with Spitzer (Sect l3.5t . 

The known contributors to the visibility measured with MIDI 
are the stellar photosphere visibility V*(/l, B), and the contribu- 
tion from the extended CSE resolved with VISIR V L (A, B). As 
this last component is already resolved spatially by a single tele- 
scope, it is also fully resolved by the interferometer, and its visi- 
bility V L (A, B) = 0. In other words, the CSE detected with VISIR 
contributes only spatially incoherent flux to the MIDI observa- 
tions. These two contributors create a visibility dependence of 
the form: 



V mo del(A, B) - 



Coherent flux density f±{A) VJA, B) 



Total flux density 



U(A) + f L (A) 



(22) 



where /*(/l) is the stellar flux density, V*(/l, B) the visibility of 
the photosphere, and /l(A) the flux density of the large CSE com- 
ponent. This can be rewritten as: 



V mo del(A, B) - 



V*(A,B) 
1 + a(A) 



(23) 



where a(A) = /z,(/l) //*(/!) is the relative contribution of the 
large component. This parameter was directly measured from 
the VISIR observations in the PAH1 and PAH2 filters (Sect.lUi. 
As we assume that ff(8 pm) = 0, this gives us the differential 
visibility expression: 



AV model (A,B) = 



AV i ,(A,B)-a(A) 
1 + a(A) 



(24) 



http://www.eso.org/sci/facilities/paranal/instruments/visir/ 



In the next paragraphs of this Section, we compare the measured 
value of AV ohs (A,B) in the PAH1 and PAH2 bands to the pre- 
dicted AV mo de\(A, B), for the two stars. 

5.3. Sensitivity domain 

Now let us discuss the sensitivity of our approach for the detec- 
tion of a compact CSE component, within the angular dimension 
and photometric contribution parameter space. We now consider, 
in addition to the large CSE component discussed precedently, a 
compact CSE contributing a fraction B(A) of the stellar flux den- 
sity, with a non-zero visibility Vc(A, B). This additional compo- 
nent represents the perturbing signal that we seek in our MIDI 
observations, once the extended VISIR CSE has been taken into 
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account. As this new component adds a partially coherent flux 
contribution, its presence will modify Eg. [231 to wards: 



Vmodel(-i. B) 



VJA,B) + B(A)V C (A,B) 
1 + a(A) + 6(A) 



(25) 



where B(A) Vc(A, B) is the coherent part of the total flux density 
of the compact CSE, and B(A) its total flux density. The sensitiv- 
ity of our study to the presence of a compact CSE can be divided 
in two domains, depending on its typical angular size pc: 

- pc < A/B: the compact component visibility Vc(A, B) will 
be significantly different from zero. In this case, the MIDI 
differential visibility measurement will show the presence 
of the compact CSE as an additional slope (to first order) 
on the differential visibility spectrum compared to the model 
expected from the VISIR extended CSE alone. However, if 
the coherent flux contribution B(A)Vc(A, B) is uniform with 
wavelength, this signature will disappear in the computation 
of the differential quantity through Eq.l2Tl 

- A/B < pc < A ID: in this case its interferometric visibility 
is negligible, and it will remain undetectable in the differen- 
tial visibility spectrum (no effect on the differential visibility 
AV m odei computed from Eq.l25ll. 

The angular resolution A/B of the interferometer in the N 
band is typically « 55 mas on the 40 m UT2-UT3 baseline, and 
ss 17 mas on the 130 m UT1-UT4 baseline, while the resolution 
of a single Unit Telescope is A/D =s 500 mas. 

It thus appears that the differential visibility technique is sen- 
sitive only to a restricted kind of compact CSE. We can only de- 
tect the color excess of the partially resolved part of a compact 
CSE component. This color excess is measured over the N band, 
relative to A = 8 pm. We are not sensitive to its absolute flux 
level. The contribution from the envelope resolved with VISIR 
is taken into account for this comparison. Within these limita- 
tions, and as shown in the following numbers, our sensitivity 
expressed as a fraction of the stellar photosphere flux is between 
1.6 and 7.5%. 



Photosphere 




L Car- UT1-UT4 



10.5 
Mum) 



11,0 11,5 12,0 12,5 13,0 13,5 




Fig. 11. Differential visibility of ^Car on the UT1-UT4 (top) and 
UT2-UT3 (bottom) baselines. The details of the different curves 
are discussed in Sect. 15.41 



5.4. i Car differential visibility spectrum 

We observed ^Car with MIDI on the UT2-UT3 baseline (B = 
40.0 m) in 2004 (#A and #C, Table H, and on the UT1-UT4 
baseline in 2006 (#E and #P). Observation #J was rejected due to 
its very low interferometric flux. The average projected baseline 
length of the two UT1-UT4 observations is B = 123.8 m. 

The differential visibility spectra are presented in Fig. QT| 
The MIDI average spectra for 2004 and 2006 are shown as thick 
black curves, with the original error bars on the absolute visibil- 
ities produced by the EWS software. The thin rising curve rep- 
resents the expected differential visibility of the photosphere of 
the star without CSE. The open squares show the expected dif- 
ferential visibility in PAH1 and PAH2 of the stellar photosphere 
plus the extended component resolved with VISIR (as defined in 
Sect.|3T2l. The width of the PAH1 and PAH2 filters bandpass are 
shown as horizontal error bars. In the upper plot of Fig. [TT] the 
contribution of the CSE to the decrease of the differential visi- 
bility in indicated for the PAH2 band. The thick grey curve is a 
linear fit to the MIDI visibility spectrum, scaled so that it equals 
zero for A = 8 pm. The thin dashed curves represent the rms dis- 
persion domain of the differential visibilities around the linear fit 
model, that equal 1.6% and 2.6% respectively for the UT1-UT4 
and UT2-UT3 baselines. 



The differential visibility of the photosphere AV*(/l, B) was 
computed from the expected limb darkened angular diameter of 
the star. For the 2004 epoch, £ Car was at its minimum diameter 
phase {(j) = 0.90), for which Kervella et al. (2004a) find a pho- 
tospheric angular diameter of #ld = 2.69 mas. The 2006 epochs 
correspond to phases of <p — 0.66 and <p = 0.46. The average an- 
gular diameter of I Car for these two phases is #ld = 3.09 mas. 
We can consider these two observations together, as the varia- 
tion of the photospheric angular diameter of the star between 
the two phases is only 66u) = +0.10 mas. This corresponds on 
the 130 m UT1-UT4 baseline to a change in visibility of 1.7% 
at 8.0yum, below the typical uncertainty on the MIDI differential 
visibilities. 

The UT2-UT3 differential visibility spectrum obtained with 
the present reduction differs in slope from the V 2 curve pre- 
sented in Fig. 6 of Kervella et al. d20061 l. The reason for this 
difference is likely to be the different definition of the detector 
masks used to extract the source spectrum and estimate the back- 
ground. This step is critical for the correct analysis of the MIDI 
data, and the computation method of the EWS software used 
for the present work (mask adapted to the observed object) ap- 
pears more robust than that of the Paris Observatory package 
(fixed mask position on the detector). But although the spectrum 
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slopes are different, the two calibrated visibility spectra agree 
well within their statistical uncertainties. From an integration of 
the measured AV(A,B), weighted by the PAH1 and PAH2 filter 
transmissions, we derive for the UT1-UT4 baseline: 



AV obs (PAHl) = -0.2 + 1 .6 %, (26) 

AV obs (PAH2) = -3.3 ± 1.6 %, (27) 
and we compute model values of: 

AV model (PAHl) = -2.8 + 0.9 %, (28) 

AV mod ei(PAH2) = -6.5 + 1.3%. (29) 
Similarly, on the UT2-UT3 baseline, we measure: 

AV obs (PAHl) = +0.1 + 2.6 %, (30) 

AV obs (PAH2) = -8.6 + 2.6 %. (31) 
and we obtain model values of: 

AV mo dei(PAHl) = -3.7 + 0.9 %, (32) 

AV mode i(PAH2) = -9.4 + 1.3 %. (33) 



As shown in Fig. [TT] and the comparison of the AV obs and 
AV m odei values above, the model differential visibilities AV mo d e i 
are smaller than the observed AV obs by w 3 ± 1.5% in PAH1 and 
PAH2. This indicates that the relative CSE flux contribution at 
8.0/im is probably of this order, as the effect of the presence of 
extended emission at this wavelength is a global negative shift of 
the AV obs curve relative to the photospheric visibility. This good 
agreement between the observed and model differential visibili- 
ties imply that the deficit observed between the expected photo- 
spheric visibility and the measured MIDI visibility is caused by 
the CSE component resolved with VISIR alone. We thus do not 
find a significant additional compact component. 

5.5. RS Pup differential visibility spectrum 

The differential visibility spectrum of RS Pup is shown in 
Fig. Q~2] This curve is the average of the three observations of 
the star labeled #H, #L and #0 in Table [2] obtained between 
phases <f> — 0.92 and <p — 0.97. The variation in projected base- 
line length is negligible between the three observations. As ob- 
served for i Car, the visibility decreases significantly from 8 to 
13//m, and the residual dispersion of the differential visibilities 
around the linear fit shown in Fig.[T2]as a thick grey line is 7.5%. 
The measured differential visibilities are well reproduced by our 
single-component extended emission model, and we obtain: 



AV obs (PAHl) = -5.6 ± 7.5 %, (34) 

AV obs (PAH2) = -14.2 ± 7.5 %, (35) 

AV mod ei(PAHl) = -2.1 ± 1.2%, (36) 

AV mod ei(PAH2) = -18.7 ± 2.0 %. (37) 



The observed and modeled values are in agreement well within 
lcr. As for £Car, we thus do not need to invoke a compact CSE 
component in addition to the VISIR resolved flux to explain the 
measured differential visibilities. Moreover, we do not detect a 
significant extended emission at 8.0 /mi at a level of 7% of the 
stellar flux, as we do not detect a global negative shift of the 
model visibilities as in the case of I Car. 
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Fig. 12. Differential visibility of RS Pup on the UT1-UT4 base- 
line. The symbols are the same as in Fig. [TT] but the vertical 
scale is different. 

5.6. N band spectroscopy 

The data acquisition of MIDI requires the observation of a vis- 
ibility calibrator that is usually chosen to be also a photometric 
standard. The absolutely calibrated spectrum of the target can be 
obtained by dividing its average MIDI spectrum by that of its 
calibrator, and then multiplying the result by the absolutely cal- 
ibrated template spectrum of the latter, for instance from Cohen 
et al. (fT999l . 

Two of the 2006 observations of i Car were obtained with 
HD 94510 as calibrator. This K1III giant is not in the Cohen 
et al. (1 19991 1 catalogue, and we derived its thermal infrared 
spectrum using its CalVir[3 limb darkened angular diameter 
(6>ld = 2.23 + 0.01 mas) and the theoretical spectrum by Castelli 
& Kurucz (2003) for an effective temperature of 5 000 K, solar 
metallicity and logg = 3.0 (parameters from Cayrel de Strobel 
et al. |2001| l. We checked a posteriori that the photometry com- 
piled by Ducati (2002) is well reproduced by the resulting model 
spectrum, in the B to L bands. This gave us the template spec- 
trum of HD94510 used for the photometric calibration of the 
I Car spectra #J and #P in Tabled 

The calibrated N band spectra of i Car is presented in Fig.[T3l 
(top). The agreement with the IRAS spectrum from Volk & 
Cohen (1989) is qualitatively satisfactory. The slightly greater 
flux in the IRAS data is due to the larger aperture over which 
the photometry has been obtained. The extended component re- 
solved with VISIR (Sect. 14.31 ) was excluded from the measure- 
ment, causing an underestimation of the total flux. There is no 
statistically significant change in the spectral energy distribution 
between the 2004 (0 = 0.90) and 2006 (0 = 0.56) epochs. The 
spectrum of RS Pup in the bottom part of Fig. [13] is the first of 
this star in this wavelength domain. It was obtained close to the 
phase of maximum light in the visible (<p = 0.9 - 1.0). Fig. [T3l 
shows that the agreement of the MIDI spectra with the broad- 
band photometry listed in Table [5] is satisfactory at the » 1 cr 
level for RS Pup (bottom curve). 

An important aspect of the MIDI spectrophotometry is that 
it is obtained with a very narrow aperture. The field of view of 
the MIDI cold stop is approximately 1.5" in radius, over which 
is superimposed a broad slit of 0.6" x 2.0", that rotates with the 



http://www.eso.org/observing/etc/ 
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Fig. 13. Thermal infrared spectrum of ^Car (top) and RSPup 
(bottom) from MIDI. The thin curve in the top figure corre- 
sponds to the 2004 observation, while the thick curve is the 
2006 data. The dashed curve is the IRAS spectrum of £ Car, and 
the open square symbols in both plots represent the photometric 
measurements listed in Table [5] excluding the extinction correc- 
tion. 



Earth rotation (due to the absence of field derotator). In the MIDI 
data processing software, the spectrum is extracted from the de- 
tector frames using a Gaussian weighting perpendicularly to the 
dispersion direction, that effectively reduces the field to « 0.4" 
FWHM. At this stage, the background is subtracted from two 
lateral bands parallel to the dispersion direction within the field 
of view, thus subtracting also the contribution of any extended 
emission located beyond a radius of 0.2 to 0.3" from the star. 
For an example of aperture photometry with MIDI on the ex- 
tended CSE of j] Car, the interested reader is refered to Chesneau 
et al. (t2005T >. 

Due to these difficulties, we did not use the MIDI spectra for 
the analysis of the CSE spectral properties in Sect. [8] and restrict 
their use to the search for spectral features that are absent from 
the N band spectra of the two stars. 



6. VINCI data 

To complete our new observations, we also consider here the ob- 
servations of ^Car obtained with the VINCI instrument in the K 
band by Kervella et al. (2006). They are well reproduced by the 
stellar photosphere and a Gaussian CSE with an angular radius 
at half maximum of 



Pi.2iim(t Car) = 2.9 + 2.3 mas 

and a relative flux contribution a of 

a 2 .2^Car) =4.2 + 0.2%. 



(38) 



(39) 



This model is of the same nature as those we adjusted to the other 
data presented here. We can therefore directly include these pa- 
rameters in our analysis. 

7. Properties of the CSEs of I Car and RS Pup 

7. 1 . Spectral energy distribution 

To retrieve the spectral energy distribution (SED) of the CSEs 
of i Car and RS Pup from the total observed flux, we need to 
estimate the contribution from the stellar photosphere. For this 
purpose, we used a combination of our new photometry in the 
near to mid-IR domains with intensity mean photometry assem- 
bled by Fouque et al. (|2007T > in the BVRJ C JHK bands from the 
literature. In the infrared, we also collected the COBE/DIRBE 
photometry of £Car at 3.5, 4.9 and 12//m (Smith et al.|2004), 
the photometry of £Car at 12, 25 and 60 fim and RS Pup at 12 
and 25 /mi from the IRAS Point Source Catalogue (IPAC [19861 ). 
We took the photometry of RS Pup at 60 and lOOyimi from the 
IRAS Small Scale Structure Catalogue (Helou & Walker [1988] ), 
as the angular size of the nebula at these wavelengths becomes 
non-negligible compared to the angular resolution of IRAS. 

The B to K photometry from Fouque et al. (2007) is al- 
ready corrected for interstellar absorption. In order to preserve 
the consistency of our data set at all wavelengths, we adopted 
their E(B - V) color excess values, namely 0.147 + 0.013 for 
( Car and 0.457 + 0.009 for RS Pup. We also used the same total- 
to-selective absorption ratio Ry = 3.23 to derive Ay extinctions 
for the two stars: A v (£Car) = 0.475 and A^RSPup) = 1.476. 
We then used these values and the interstellar extinction curves 
in the mid-IR by McClure (2008) to correct all our photomet- 
ric measurements between 3 and 30 fim. As a remark, the curves 
from this author are in excellent agreement with the independent 
results by Chapman et al. ( 120081 ), so we do not expect any sig- 
nificant systematics from this operation. Longwards of 30/im, 
the correction to apply to the flux density is negligible (below 
1 % for t Car and 2% for RS Pup), compared to the uncertainties. 
The complete set of measurements is summarized in Table [5] 

To model the photospheric emission, we used tabulated stel- 
lar atmosphere models obtained with the ATLAS9 simulation 
code that we retrieved from Castelli & Kurucz's (2003 ) library. 
The grid we chose was computed for solar metallicity and a tur- 
bulence velocity of 2km/s. We finely interpolated this grid in 
order to compute spectra for any effective temperature and any 
surface gravity. We multiplied the spectrum by the solid angle 
of the stellar phosphere, n6^ D /4, where #ld is the limb darkened 
angular diameter, in radians, in order to scale the model. The 
photometric data were adjusted to the model taking into account 
the spectral response of each instrument listed in Table [5] Since 
we assume that there is no detectable excess (5% or less) below 
the 2.2 fim K band, we included all the photometric measure- 
ments bluer than the K band to estimate the angular diameter and 
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the effective temperature. Note that we did not adjust the surface 
gravity, since the broadband photometry is mostly insensitive to 
this parameter. 

Although we are not interested here in an unbiased estimate 
of the angular diameter #ld, the best-fit values of 6*LD(^Car) = 
2.92 mas is only 2% smaller than the mean angular diameter of 
6 LD (£ Car) = 2.99+0.01 mas measured by Kervellaetal. (12004d l 
in the K band. For RS Pup, the agreement between the present 
best-fit value of 6>ld(RS Pup) = 0.93 mas and the 1.02 mas pre- 
diction by Moskalik & Gorynya (2006) is also reasonably good. 
The best-fit effective temperatures T e g of £Cai and RS Pup are 
respectively 4860 K and 5040 K. One should note however that 
there is a correlation between T e ff and the angular diameter in 
this fit. Following Kervella et al. (2004c), we consider an effec- 
tive gravity of log g = 1.5, typical of a long period Cepheid, but 
changing this parameter by + 1.0dex does not significantly affect 
the quality of the fit, nor the amount of infrared excess. 

The fitting process produces a good agreement between the 
model spectrum and the BtoK band photometric measurements, 
with a reduced^- 2 of 0.7 for t Car, and 2.2 for RS Pup. It should 
also be noted that in the case of RS Pup, the presence of the neb- 
ulosity makes the estimation and correction of the interstellar 
and circumstellar reddening quite difficult, which might explain 
in part the slightly worse agreement. The result of this fit is pre- 
sented in Fig. [14]for £ Car and RS Pup. In both cases, the pres- 
ence of excess radiation in the thermal infrared domain is clearly 
detected. 



7.2. Spectral and spatial structure of the infrared excess 

We discuss in this paragraph the distribution of the infrared ex- 
cess as a function of wavelength, both in the spectral energy 
distribution and as spatially extended emission. The different 
properties of the CSEs of ^Car and RS Pup are summarized in 
Fig. Q3] While we discuss here the overall properties of the ob- 
served circumstellar emission, we postpone the presentation of a 
detailed model of the Cepheid CSEs to a future publication, that 
will incorporate the present results and the scattered light obser- 
vations of RS Pup's nebula obtained by Kervella et al. (2008). 

To convert the angular values determined previously into lin- 
ear distances from the star, £Cai, we used the distance of 566 pc 
measured interferometrically by Kervella et al. (2004a). While 
we do not resolve the CSE in the Spitzer/MlPS images at 24 and 
70 /mi, the angular resolution of these images allows us to set 
an upper limit on its extension of « 10" and 20", respectively. 
The shortest wavelength at which the excess emission of £ Car 
is detected is « 8.0//m. For RSPup, we assume a distance of 
1 992 pc (Kervella et al. 120081 . We included in Fig.[l5]panel (f) 
the extension of the RS Pup nebula measured with IRAS at 60 
and 100 /mi, respectively of 2.1 and 3.7 arcmin in angular diam- 
eter (DeasyQIHl. 

The CSEs of ^Car and RS Pup appear different in many re- 
spects. The distance corrected flux density of £ Car's envelope 
appears as continuously decreasing with wavelength ( Fig. [15] 
panel (c)), with an almost exponential dependence. The obser- 
vation of a CSE flux contribution at 2.2 //m by interferometry 
(Kervella et al. 2006) shows that a very hot and compact com- 
ponent is present in its CSE. From simple radiative equilibrium 
considerations and its measured linear extension (Fig. [15] panel 
(e)), this implies that the envelope starts within a few stellar radii, 
and extends continuously up to at least several hundred astro- 
nomical units from the star. 

The envelope of RS Pup is on average much larger and colder 
than that of I Car. The most striking difference between £ Car and 
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Fig. 14. Synthetic spectra of I Car and RS Pup (solid line) with 
the photometric measurements listed in Table[5](dots). The resid- 
uals of the subtraction of the synthetic spectrum from the mea- 
surements is shown in panels (c) and (d) of Fig.fTBI 



RS Pup's CSEs is their brightness around 70/mi (Fig. [15] panels 
(c) and (d)). While the two stars exhibit approximately the same 
brightness around 10/mi, RS Pup's envelope is « 2 000 times 
brighter than £ Car's at 70 /mi. The presence of such a consider- 
able emission at 70 //m indicates that the temperature of the ma- 
jority of the dusty material in the CSE of RS Pup is around 40 K, 
as already established by McAlary & Welch (1986). However, 
our VISIR and MIDI observations around 10 /mi as well as the 
IRAS photometry at 12 /mi show that a warm component is also 
present in the CSE of RS Pup. The angular extension of this 
warm component, as well as its flux contribution, appear to be of 
the same order of magnitude as that of ^Car, as shown in panels 
(c) to (f) of Fig. [15] The observed properties of RS Pup's CSE 
therefore suggest that it consists of two components: 

- A warm and compact component, located at typical distances 
of a few 100 to a few 1 000 AU from the star, with a typical 
temperature of a few 100 K (maximum flux around 10 /mi), 

- A cold component distributed within a relatively large hol- 
low ellipsoid with a typical radius of a few 100000 AU and 
temperature around 40 K (maximum flux around 70/mi). 
This cold component is also the scattering medium for the 
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Fig. 15. Comparison of the spatial and spectral properties of the CSEs of ^Car (left) and RSPup (right). In all panels, the open 
squares represent the CSE resolved with VISIR and the open circles the CSE of IQw in the K band (resolved with VINCI), (a) 
and (b): Excess flux density relatively to the photospheric emission, with a different vertical scale for t Car and RS Pup. (c) and 
(d): Flux density from the CSE alone (stellar contribution subtracted), multiplied by the squared distance of each star, (e) and (f): 
Half-maximum linear radius of the envelope as a function of wavelength. The photospheric linear size is shown as a solid line. The 
Spitzer observations of I Car at 24 and 70 fim provide upper limits only, as the CSE was not resolved. Ad hoc models are shown as 
thin dashed curves to guide the eye on the different plots. 
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light echoes observed by Havlen ( 119721 1 and Kervella et 
al. ([2008). It could also host the thermal infrared "heat 
echoes" predicted by Mayes et al. (1985 ), although their ob- 
servation appears extremely challenging. 

In contrast to RS Pup, the CSE of £Car appears to be devoid of 
a cold dust component, but its warm component (resolved with 
VISIR and MIDI) is comparable to RS Pup's, with a typical size 
of 500 AU and a temperature of a few 100K. 

8. Discussion 

We present in this Section our interpretation of the properties of 
the CSEs of ^Car and RS Pup. We then examine the mass loss 
scenarios predicted by existing Cepheid models, in view of these 
properties. 

8. 1 . Origins and evolution of the CSEs 

In spite of their relatively similar physical properties, £ Car and 
RS Pup appear to be surrounded by two rather different kinds 
of CSEs. The warm (lOyum) components of the CSEs of both 
stars are probably created by a continuous mass loss from the 
Cepheids themselves, likely driven by pulsation as discussed by 
Merand et al. (120071) . The hot CSEs observed in the K band 
around £Car (Sect. [6), around Polaris and 5 Cep (Merand et 
al. I20061 . and around Y Oph (Merand et al. 120071 ) very likely 
represent the innermost emission (within a few stellar radii of 
the Cepheid) of the warm component we observed with VISIR 
and MIDI around our two targets. Along the same direction, 
Bohm-Vitense & Love ( 1994) suggested that the fixed (non pul- 
sating) Ha absorption feature that they observe in the spectrum 
of t Car is caused by a CSE with a size of the order of 1000 AU. 
This dimension is in good agreement with the typical size of 
the compact component we detected. More recently, Nardetto et 
al. (12008 ) confirmed the presence of a non-pulsating component 
in the H a profiles of both £ Car and RS Pup. 

However, the large ellipsoidal cold dust shell surrounding 
RS Pup is probably made of compressed interstellar material. 
We speculate that the mechanism for this compression could 
be linked to the past evolution of RS Pup. Its evolutionary mass 
was determined by Caputo et al. d20051 > to be 12.8 + 0.9 M , 
and it thus spent its lifetime on the main sequence as a B2V 
hot dwarf, most probably fast rotating, with an effective tem- 
perature around 20 000 K. A good example of such a Cepheid 
precursor is the nearby Be star Achernar. This spectacularly flat- 
tened fast rotating B3 dwarf was shown recently to experience 
mass loss through a fast stellar wind ejected from its overheated 
polar caps (Kervella & Domiciano de Souza 2006). Such a pref- 
erential mass loss along the polar axis of the star could explain 
the elongation observed in the 70 pm image of RS Pup presented 
in Fig. [3] The interstellar medium (ISM) would be more effi- 
ciently compressed along the polar axis, thus elongating the cav- 
ity carved by the star in the dusty ISM. Another additional pres- 
sure would come from the enhanced ultraviolet flux as seen from 
above the stellar poles, compared to the darkened equatorial belt 
(due to the Von Zeipel effect: see Aufdenberg et al. 2006 for an 
application to Vega). 

A plausible analog of a "young" RS Pup is provided by the 
Herbig Be object HD 200775 and its surrounding dusty reflec- 
tion nebula VdB 139 (also referred to as NGC 7023). The cen- 
tral object of this nebula is a binary with a composite spec- 
trum B2Ve. The primary component has a mass of 1 1 M and 
the secondary component 9M , with a system age of only 



~ lOOOOOyears (Alecian et al. |20081 >. Considering their masses, 
the two stars will cross the instability strip in the HR diagram and 
become Cepheids within a few tens of million years. The bright 
pair illuminates a large circumstellar nebula that was extensively 
studied at visi ble (Wi tt et al. 120061 Berne et al. 120081 1. infrared 
(Werner et al. I2004bl Sell gren e t al. 120071 1. and submillimetric 
wavelengths (Gerin et al. [19981 . The distance of HD 200775 
is uncertain. It was estimated at 430 pc by van den Ancker et 
al. ([T997 ) based on the Hipparcos parallax, before the star was 
recognized as a binary. Based on Wolf diagrams, Kun (1998 ) 
find a closer distance of 300 pc. We adopt here the value from 
van den Ancker et al. ([1997 ), although it may be biased by un- 
recognized orbital motion. The angular radius of the reflection 
nebula of about 10 arcminutes correponds to a linear extension 
of the order of 200 000 AU from the central star, which is com- 
parable to the linear radius of the RS Pup nebula. In spite of its 
young age, HD 200775 already started to carve a cavity in the 
interstellar dust, and this cavity presents a biconical structure 
that is well defined at visible, infrared (Werner et al. 2004bJ and 
radio wavelengths (Fuente et al. 1998). A visible image from 
the red POSS-II atlas and Spitzer/MlPS images of VdB 139 at 
24 and 70 pm (retrieved from the SSC archive) are shown in 
Fig-US] to be compared with the images of RS Pup presented in 
Fig. [3] The overall distribution of the thermal infrared emission 
in VdB 139 appears similar to RSPup, with irregular filamen- 
tary structures and knots. These structures, that are also observ- 
able in the visible images of the nebula in scattered light (Witt et 
al. 2006 ), are typical of the turbulent ISM in star forming regions 
(Boulanger[3nn71 Dwarkadas|2008]!. 

8.2. Evolutionary vs. pulsational masses of Cepheids 

From a model of the dust emission based on IRAS photometry, 
Havlen d 19721 ) and McAlary & Welch i 19861 1 obtain a dust mass 
of the nebula surrounding RSPup of approximately 0.03 M , 
converting to a total gas+dust mass of the order of 3M Q . If the 
cold component of RS Pup's nebula is indeed interstellar in ori- 
gin, the actual mass originating from the star itself would be 
much smaller. Such a scenario would strengthen the theoretical 
results by Caputo et al. (120051 1 and Neilson & Lester ([2008a 
2008b) that predict a higher mass loss rate for short period 
Cepheids. On the other hand, if the totality of the mass of the 
RS Pup nebula originates from the star, this would be in con- 
tradiction with the modeling by Caputo et al. (2005 ), who es- 
tablished that the pulsation and evolutionary masses of RS Pup 
agree well, without evolutionary mass loss. 

The conclusion of Caputo et al. (2003J is in striking con- 
tradiction with the computations by Keller (2008) using the 
Padova evolutionary modeling code. This author predicts a uni- 
form mass discrepancy between the evolutionary and pulsation 
masses, the former being 17 + 5% larger than the latter, inde- 
pendently of the pulsation period. Applied to RS Pup's mass 
of « 13 M , this converts into a total mass loss of 2.2 M Q , in 
excellent agreement with the 3 M nebular mass computed by 
McAlary & Welch ( TT9861 . Keller ( 120081 a ttribut es the discrep- 
ancy in derived mass with Caputo et al. (2005) to a different 
treatment of convective overshoot between the two models. 

However, a problem with Keller's (2008 ) model is that the 
application of the same 17% mass loss fraction to £ Car's (whose 
mass is also ss 13M ) leads naturally to the same « 2M total 
mass loss, and consequently =s 0.02 M Q of dust. Such a consid- 
erable mass should be observable as a CSE similar to RS Pup 
at 24 and 70 fim. According to the period-age relationship es- 
tablished by Bono et al. (2005), I Car and RS Pup are between 
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Fig. 16. Structure of VdB 139 at visible and infrared wavelengths: POSS-II red image (left), Spitzer/MTPS images at 24 pm (center) 
and at 10 pm (right). The position of the star is marked with and the linear scale of the image is given for a distance of 430 pc. 



17 and 19 million years old. As the main sequence lifetime of a 
13 M star is « 15 to 17Myr, and assuming that the mass is lost 
mostly during the He burning phase, the star would have to ex- 
pel of the order of 10~ 6 M Q /year continuously over 2 Myr. Such a 
rapid and extremely high mass loss rate would necessarily create 
a massive nebula around /'Car, that is not observed at present. A 
differential dissipation of their CSEs in the ISM is also difficult 
to invoke based on the short timescales involved. 



9. Conclusion 

Thanks to a combination of instrumentation in the thermal in- 
frared covering a range of spatial resolutions from arcminutes to 
milliarcseconds, we presented a characterization of the CSEs of 
I Car and RS Pup both spatially and spectrally. The detected in- 
frared excess relative to the stellar photospheric emission ranges 
between 10% and 50% for /'Car, from 8 to 70 /an. A similar 
level of emission is also observed around RS Pup, and a much 
more considerable contribution fully outshines the stellar emis- 
sion longwards of 24 /an. We thus conclude that both CSEs host 
a warm component with a similar brightness and spatial exten- 
sion around 10 /mi, and that RS Pup contains in addition a cold 
dusty ellipsoidal shell located at a few 100 000 AU from the star. 
This cold component is not detected around (Car. 

Our interpretation is that the warm CSE components are the 
consequence of ongoing mass loss from the stars. Based on its 
spectral and spatial properties, the cold dust shell of the CSE of 
RS Pup is most likely of interstellar origin. We speculate that it 
was shaped in its present elongated ellipsoidal form during the 
main sequence lifetime of the then fast rotating B2V progeni- 
tor of RS Pup. The ISM was compressed by the strong stellar 
wind and UV radiation preferentially along the polar axis of the 
star. We also propose that the reflection nebula surrounding the 
Herbig Be star HD 200775 is an analogue of RS Pup shortly after 
it formed as a star. In this scenario, the absence of a cold dusty 
shell around i Car would be explained by a lower ISM density 
around this star at the time of its formation. Conversely, the pres- 
ence of very large dusty envelopes such as that of RS Pup would 
require a high density of dust in the ISM at the time of the forma- 
tion of the Cepheid progenitor. In this context, their occurrence 
around Cepheids is probably relatively rare. 

The presence of warm CSEs around the two brightest long- 
period Cepheids shows that many long-period members of this 
class, if not all, could be surrounded by warm circumstellar en- 



velopes created by ongoing mass loss. This is also the conclu- 
sion of Merand et al. (2007) from interferometric observations in 
the near-IR. In the presently studied cases of I Car and RS Pup, 
the warm CSE flux contribution is negligible in the visible and 
limited in the near-IR (< 5% of the photosphere). Their direct 
photometric impact on the Cepheid distance scale at these wave- 
lengths will thus be minor or negligible. However, as shown 
by Neilson et al. (2008), the differential behavior of long and 
short period Cepheids with respect to mass-loss could affect the 
slope of the P-L relations. Although weak in the near-infrared, 
CSEs also have an influence on the determination of Cepheid 
distances using the Baade-Wesselink method, and particularly 
its interferometric version (see e.g. Merand et al. 120061 1 . If not 
properly taken into account, the spatially incoherent CSE flux 
contribution biases the interferometric measurement of the an- 
gular diameter, and the derived distance as a consequence. The 
detected mid-IR excess is more considerable (x: 15% at 10pm) 
and could noticeably complicate the definition of thermal in- 
frared P-L relations. As suggested by Ngeow & Kanbur (2008), 
such relations nevertheless appear highly desirable for future ex- 
tragalactic Cepheid observations with e.g. the James Webb Space 
Telescope. 
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Table 5. Overview of the infrared photometry of I Car and RS Pup from the present work and bibliography. All flux densities are 
corrected for interstellar extinction, except longwards of 30/~m, where it is considered negligible. The "•" symbols in the "Fit" 
column indicate the photometry used for the synthetic spectrum fitting, /lo is the average of the half maximum extreme wavelengths 
of each instrument/filter combination, and AA the half-bandwidth. A e ff is the effective wavelength of the measurement considering 
the spectral energy distribution of the star. The "Ap." column lists the effective radius of the aperture over which the flux density was 
measured, when this parameter was available. In the "Ref." column, the references are the following: F07 — Fouque et al. (2007), 
186 = IPAC (fT986l H88 = Helou & Walker ([19881 , K08 = present work, S04 = Smith et al. ( 12004. . "a" is the measured flux density 
excess, expressed in percentage of the photospheric flux (see Sect. 17. 1 I for details), except for RS Pup's three longest wavelengths, 
where it is expressed in stellar flux units, followed by "x". The last column "Ncr" gives the detection level in number of times the 
RMS uncertainty. 
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18.8 ±5.2 


3.6 


IRAS** 


25 fim 


24 ±6 


23.34 


25 


3.96 ± 0.28 


10" 


15 


0.72 ± 0.05 


186 


96 ± 14 


7.0 


MIPS 


Ml 


23.5 ±2.7 


23.68 


180 


4.21 ±0.17 


10" 


15 


0.79 ± 0.03 


K08 


145 ± 9.8 


15 


IRAS 


60 jjm 


62 ± 17 


59.32 


CO 


1.19 ±0.26 


10" 


14 


13.9 ±3.1 


H88 


220 ± 49x 


4.5 


MIPS*** 


M2 


70.4 ± 9.5 


71.42 


oo 


9.94 ± 0.50 


10" 


15 


16.9 ±0.9 


K08 


463 ± 23x 


20 


IRAS 


100/im 


103 ± 18 


100.3 


oo 


7.87 ± 1.57 


io- 


15 


26.4 ±5.3 


H88 


1 430 ± 290x 


5 



* The angular resolution of COBE/DIRBE is w 0.7°, its aperture is thus considered as infinitely large. 

** The IRAS Point Source Catalogue (IPAC 1986 1 apertures were taken as the 80% encircled energy from Beichman et al. ( 1988 1. 
*** The MIPS/M2 flux density of RS Pup's CSE is the integral of the Gaussian profile adjusted in Sect. 13.31 



